Entanglement distribution over the subsystems and its invariance 
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We study the entanglement dynamics of two qubits, each of which is embedded into its local 
amplitude-damping reservoir, and the entanglement distribution among all the bipartite subsystems 
including qubit-qubit, qubit-reservoir, and reservoir-reservoir. It is found that the entanglement can 
be stably distributed among all components, which is much different to the result obtained under 
the Born-Markovian approximation by C. E. Lopez et al [Phys. Rev. Lett. 101, 080503 (2008)], 
and particularly it also satisfies an identity. Our unified treatment includes the previous results as 
special cases. The result may give help to understand the physical nature of entanglement under 
decoherence. 

PACS numbers: 03.65.Ud, 03.65.Yz, 42.50.Lc 



I. INTRODUCTION 

Entanglement in quantum multipartite systems is a 
unique property in quantum world. It plays an important 
role in quantum information processing [l|, l2|- Therefore, 
the study of its essential features and dynamical behav- 
ior under the ubiquitous decoherence of relevant quan- 
tum system has attracted much attention in recent years 
|3l-[l8|. For example, it was found that the entanglement 
of qubits under the Markovian decoherence can be termi- 
nated in a finite time despite the coherence of single qubit 
losing in an asymptotical manner [3]. The phenomenon 
called as entanglement sudden death (ESD) [4, 5] has 
been observed experimentally p, 0] . This is detrimental 
to the practical realization of quantum information pro- 
cessing using entanglement. Surprisingly, some further 
studies indicated that ESD is not always the eventual fate 
of the qubit entanglement. It was found that the entan- 
glement can revive again after some time of ESD l8^llQ| , 
which has been observed in optical system [ll|, [l2|- It 
has been proven that this revived entanglement pla ys a 
constructive role in quantum information protocols [10|. 
Even in some occasions, ESD does not happen at all, in- 
stead finite residual entanglement can be preserved in the 
long time limit |13l-[l6| . This can be due to the structured 
environment and physically it results from the formation 
of a bound state between the qubit and its amplitude 
damping reservoir [13, [l8| . These results show rich dy- 
namical behaviors of the entanglement and its characters 
actually have not been clearly identified. 

Recently, Lopez et al. asked a question about where 
the lost entanglement of the qubits goes [i9|- Inter- 
estingly, they found that the lost entanglement of the 
qubits is exclusively transferred to the reservoirs under 
the Markovian amplitude-damping decoherence dynam- 



* Electronic address: |anjhong@lzu.edu.cn| 
"I" Electronic address: phyohch@nus.edu. sg' 



ics and ESD of the qubits is always accompanied with 
the entanglement sudden birth (ESB) of the reservoirs. 
A similar situation happens for the spin entanglement 
when the spin degree of freedom for one of the two parti- 
cles interacts with its momentum degree of freedom [20| . 
All these results mean that the entanglement does not 
go away, it is still there but just changes the location. 
This is reminiscent of the work of Yonac et al. [21], in 
which the entanglement dynamics has been studied in 
a double Jaynes-Cummings (J-C) model. They found 
that the entanglement is transferred periodically among 
all the bipartite partitions of the whole system but an 
identity (see below) has been satisfied at any time. This 
may be not surprising since the double J-C model has no 
decoherence and any initial information can be preserved 
in the time evolution. However, it would be surprising 
if the identity is still valid in the presence of the deco- 
herence, in which a non-equilibrium relaxation process is 
involved. In this paper, we show that it is indeed true for 
such a system consisted of two qubits locally interacting 
with two amplitude-damping reservoirs. It is noted that 
although the infinite degrees of freedom of the reserviors 
introduce the irreversibility to the subsystems, this result 
is still reasonable based on the fact that the global sys- 
tem evolves in a unitary way. Furthermore, we find that 
the distribution of the entanglement among the bipartite 
subsystems is dependent of the expHcit property of the 
reservoir and its coupling to the qubit. The rich dynam- 
ical behaviors obtained previously in the literature can 
be regarded as the special cases of our present result or 
Markovian approximation. Particularly, we find that, in- 
stead of entirely transferred to the reservoirs, the entan- 
glement can be stably distributed among all the bipartite 
subsystems if the qubit and its reservoir can form a bound 
state and the non-Mar kovian effect is important, and the 
ESD of the qubits is not always accompanied with the 
occurrence of ESB of reservoirs. Irrespective of how the 
entanglement distributes, it is found that the identity 
about the entanglement in the whole system can be sat- 
isfied at any time, which reveals the profound physics of 



the entanglement dynamics under decoherence. 

This paper is organized as follows. In Sec. |TT1 the 
model of two independent qubits in two local reservoirs is 
given. And the dynamical entanglement invariance is ob- 
tained based on the exact solution of the non-Markovian 
decoherence dynamics of the qubit system. In Sec. HIH 
the entanglement distribution over the subsystems when 
the reservoirs are PBG mediums is studied explicitly. A 
stable entanglement-distribution configuration is found 
in the non-Markovian dynamics. Finally, a brief discus- 
sion and summary are given in Sec. IIVI 



II. THE MODEL AND THE DYNAMICAL 
ENTANGLEMENT INVARIANCE 

We consider two qubits interacting with two uncorre- 
lated vacuum reservoirs. Due to the dynamical indepen- 
dence between the two local subsystems, we can firstly 
solve the single subsystem, then apply the result obtained 
to the double-qubit case. The Hamiltonian of each local 
subsystem is [22] 
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where cr± and ujq are the inversion operators and transi- 
tion frequency of the qubit, a\ and ak are the creation 
and annihilation operators of the /c-th mode with fre- 
quency ujk of the radiation field. The coupling strength 
between the qubit and the reservoir is denoted by gk = 
cjoe/c -d/v^Sio^jfel/, where e/c and V are the unit polariza- 
tion vector and the normalization volume of the radiation 
field, d is the dipole moment of the qubit, and Sq is the 
free space permittivity. 

For such a system, if the qubit is in its ground state |— ) 
and the reservoir is in vacuum state at the initial time, 
then the system does not evolve to other states. When 
the qubit is in its excited state |+), the system evolves 
as 

m)) = bit) \+, {o}fc) + J2 hit) I-, {i}fe) • (2) 



Here | — , {l}/c) denotes that the qubit jumps to its ground 
state and one photon is excited in the k-th mode of the 
reservoir. b{t) satisfies an integro-differential equation 



b{t) + iuJob{t) 
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where the kernel function f{t — r) = J^ d(jjJ{(jj)e '^^^^ ^^ 
is dependent of the spectral density J(cj) = 
^k\9k\ ^(^ — ^k)- Introducing the normalized col- 
lective state of the reservoir with one excitation as 
|i) = TT-T X^/j. ^/c(^) |{l}fe) and with zero excitation as 

\^)r = l{0}fc) [^' E^- c^^ b^ written as \(j){t)) = 

h(t) 1+) |6)^ + h{t) I-) |i)^, where h{t) = Jl-\h{t)\\ It 



should be emphasized that the introducing of normal- 
ized collective state is not a reduction of present model 
to the J-C model [21|, as noted in [l9|. The dynamics is 



given by Eq. (jSj, which is difficult to obtain analytically 
since its non-Markovian nature. In general the numerical 
integration should be used. 

It is emphasized that our treatment to the dynamics of 
the system is exact without resorting to the widely used 
Born-Markovian approximation. To compare with the 
conventional approximate result, we may derive straight- 
forwardly the master equation from Eq. (|2]) after tracing 
over the degree of freedom of the reservoir [23| , 

p{t) = -zA(t)[a+a_,p(t)]+r(t)[2a_p(t)a+ 

-a+a_p(t)-p(t)a+a_], (4) 



where the time-dependent parameters are given by 
A{t) = -Im[^], r{t) = -Re[^]. The time-dependent 
parameters A{t) and T{t) play the roles of Lamb shifted 
frequency and decay rate of the qubit, respectively. The 
integro-differential equation ([3j) contains the memory ef- 
fect of the reservoir registered in the time- nonlocal kernel 
function and thus the dynamics of qubit displays non- 
Markovian effect. If the time- nonlocal kernel function 
is replaced by a time-local one, then Eq. dH recovers 
the conventional master equation under Born-Markovian 
approximation [24]. 

According to the above results, the time evolution of a 
system consisted of two such subsystems with the initial 
state 1^(0)) = {a |-, -) + /3 |+, +)) |0> |6> is given by 



|^(t))=a|-,6)J-,6)2+/3|(/>(t))J(/>(t)), 



(5) 



where a and P are the coefficients to determine the initial 
entanglement in the system. From p = |<l>(t))(<I>(t)|, one 
can obtain the time-dependent reduced density matrix of 
the bipartite subsystem qubit l-qubit2 {qiq2) by tracing 
over the reservoir variables. It reads 
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where p= |/36(t)|^ 6(t)2 and x = 1 - |/3p|6(t)|'^ - 2p. Sim- 
ilarly, one can obtain the corresponding reduced density 
matrices for other subsystems like reservoir 1-reservoir 2 
(rir2) and qubit-reservoir (qiTi, qir2, <72^i, <72^2). 

Using the concurrence [25] to quantify entanglement, 
we can calculate the entanglement of each subsystem as 
Cm = max{0,(5m} with Qm for different bipartite parti- 
tions labeled by m as 

Q,,,,=2\al3\\b{t)f-2p, (7) 

Qr,r,=2\af3\b{tf-2p, (8) 

Q,,r^ = m^mm = Qq2r,. (9) 

Qq,r, = 2\a/3b{t)\b{t) -2p = Qq,r,- (10) 



One can verify that Qm in Eqs. (|7))- (fTU)) satisfy an iden- 
tity 
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2Q,,r,=2\aP\, (11) 



where 2|a/3| is just the initial entanglement present in 
qiq2- Eq. (pTj) recovers the explicit form derived in a 
double J-C model [21] when each of the reservoirs con- 
tains only one mode, i.e. J{uj) = g'^6{uj — cjo), where the 
decoherence is absent and the dynamics is reversible. It 
is interesting that this identity is still valid in the present 
model because the reservoirs containing infinite degrees 
of freedom here lead to a completely out-of-phase inter- 
action with qubit and an irreversibility. Furthermore, 
one notes that the identity is not dependent of any detail 
about 6(t), which only determines the detailed dynami- 
cal behavior of each components in Eq. (pT]) . This result 
manifests certain kind of invariant nature of the entan- 
glement. 

Eq. (pT]) can be intuitively understood by the global 
multipartite entanglement of the whole system. The 
global entanglement carried by the subsystem (qiVi) 
(^2^2) can be straightforwardly calculated from Eq. (|5]) 
by generalized concurrence [26] as 2|Qf/3|, which, coincid- 
ing with the bipartite entanglement initially present in 
qiq2, just is the right hand side of Eq. (fTTj) . Since there 
is no direct interaction between (qiVi) and (^2^2), this 
global entanglement is conserved during the time evo- 
lution. From this point, ojir result is consistent with 
the one in Refs. \]3\ and [21]. Another observation of 
Eq. ([TT]) is that the different coefficients in the left hand 
side are essentially determined by the energy/information 
transfer among the local subsystems. Explicitly, in our 
model the total excitation number is conserved, so the 
energy degradation in qi with factor b{t) is compensated 
by the energy enhancement in r^ with factor b{t). This 
causes that Qq^q^-, Qrir2-> ^ud Qq^r2i i^ all of which the 
double excitation is involved, have similar form except 
for the different combinations of h(t) and h{t) in Eqs. 
(O, ([8]), and (p^Qj) . The dynamical consequence of the 
competition of the two terms in these equations causes 
the sudden death/birth of entanglement characterized by 
the presence of negative Q. A different case happens for 
Qqiri^ where only single excitation is involved and no 
sudden death is present. With these observation, one can 
roughly understand why such combination in left hand 
side of Eq. (fTTj) gives the global entanglement. 

The significance of Eq. (pT]) is that it gives us a guide- 
line to judge how the entanglement spreads out over all 
the bipartite partitions. It implies that entanglement 
is not destroyed but re-distributed among all the bipar- 
tite subsystems and this re-distribution behavior is not 
irregular but in certain kind of invariant manner. The 
similar invariant property of entanglement evolution has 
also been studied in Ref. [275. 

In the following we explicitly discuss the entanglement 
distribution, especially in the steady state, by taking the 
reservoir as a photonic band gap (PBG) medium [28. 



CqiqiOA 




FIG. 1: Entanglement evolutions of each bipartite partitions 
for the case of cjq < cJc- The parameters used are cjq = O.lcJc 
and r] = 0.2. 



|29J and compare it with the previous results. We will 
pay our attention mainly on the consequence of the non- 
Markovian effect on the entanglement distribution and 
its differences to the results in Refs. [19.] and |21| . 



III. ENTANGLEMENT DISTRIBUTION IN PBG 
RESERVOIRS 

For the PBG medium, the dispersion relation near the 
upper band-edge is given by ^30|] 



Uk = Ldc ^ A{k - /co)^. 



(12) 



where A ^ ujc/k^^ ujc is the upper band-edge frequency 
and ko is the corresponding characteristic wave vector. 
In this case, the kernel function has the form 



fit-T)=r] I —e-'^'^'-^Uk, 



^k 



(13) 



where rj = 



ou^d^ 



is a dimensionless constant. In solv- 



ing Eq. (|3]) for 6(t), Eq. ([13]) is evaluated numerically. 
Here we do not assume that k is replaced by ko outside 
of the exponential |13Ml5l l3l| . So our result is numer- 
ically exact. In the following we take cUc as the unit of 
frequency. 

In Figs. [T]and[2l we show the entanglement evolutions 
of each subsystem for two typical cases of ujq < uOc and 
uoq > uOc^ which correspond to the atomic frequency be- 
ing located at the band gap and at the upper band of 
the PBG medium, respectively. In the both cases the 
initial entanglement in qiq2 begins to transfer to other 
bipartite partitions with time but their explicit evolu- 
tions, in particular the long time behaviors, are quite 
different. In the former case, the entanglement could 
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ulation in the long-time limit. This phenomenon known 
as population trapping [34] is responsible for the sup- 
pression of the spontaneous emission of two-level system 
in PBG reservoir and has been experimentally observed 
[28|, [29I, [35| . Such population trapping just manifests the 
formation of bound states between qi and r^ [18], which 
has been experimentally verified in [35]. Consequently, 
qi and r^ are so correlated in the bound states that the 
initial entanglement in qiq2 cannot be fully transferred 
to rir2. The oscillation during the evolution is just the 
manifestation of the strong non-Markovian effect induced 
by the reservoirs. On the contrary, if uoq is located in the 
upper band, then |6(oo)p = and the qubits decay com- 
pletely to their ground states. In this case the bound 
states between qi and r^ are absent and, according to 
Eq. (pTj). the initial entanglement in qiq2 is completely 
transferred to rir2, as clearly shown in Eq. (J8j). 



FIG. 2: Entanglement evolutions of each bipartite partitions 
for the case of cjq > ojc- The parameters used are cjq = lO.OcJc 
and Tj = 0.2. 



be distributed stably among all possible bipartite parti- 
tions. Fig. [TJa) shows that after some oscillations, a size- 
able entanglement of qiq2 is preserved for the parameter 
regime of 0.3 < a < 1. Remarkably, the entanglement 
in qiViii = 1,2) forms quickly in the full range of a [Fig. 
Wis^)] and dominates the distribution. On the contrary, 
only slight entanglement of rir2 is formed in a very nar- 
row parameter regime 0.6 < a < 1, as shown in Fig. 
[TJb). However, when uoq is located at the upper band 
of the PBG medium, the initial entanglement in qiq2 is 
transferred completely to the rir2 in the long-time limit, 
as shown in Fig. [2l At the initial stage, qiViii = 1, 2) and 
^1^2 (^2^1) are entangled transiently, but there is no sta- 
ble entanglement distribution. This result is consistent 
with that in Refs. ^, [32|, [33|] . It is noted that the entan- 
glement in qiVi comes from two parts: one is transferred 
from gig2, the other is created by the direct interaction 
between qi and r^. This can be seen clearly from Fig. [TJc) 
when a is very small, the initial entanglement of qiq2 is 
very small, while that of qiVi is rather large, which just 
results from the interaction between qi and r^. Another 
interesting point is a stable entanglement can even be 
formed for the non- interacting bipartite system qir2 [see 
[U^d) when 0.5 ^ ol < 1]. This entanglement transfer also 
results from the local interaction between qi and r^ flQj . 

It is not difficult to understand these rich behaviors of 
entanglement distribution according to Eqs. ([7|)- (p!Q|) and 
its invar iance (fTTj) . From these equations, one can clearly 
see that the entanglement dynamics and its distributions 
in the bipartite partitions are completely determined by 
the time-dependent factor \h(t)\^ of single-qubit excited- 
state population. Fig. [3] shows its time evolutions for the 
corresponding parameter regimes presented above. We 
notice that \h{oo)\^ 7^ when ujq is located at the band 
gap, which means that there is some excited-state pop- 
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FIG. 3: Time evolution of time-dependent factor of the 
excited-state population for two parameter regimes ujq = 
O.lcJc (solid line) and lO.OcJc (dashed line). 77 is taken as 0.2. 
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FIG. 4: (Color online) Entanglement evolution when a — 
\l\f2 (a), Of = 0.57 (b), and a = 0.28 (c). The parameters 
used here are the same as Fig. [1] 

In addition, in Refs. [19|, [32|, l33| it was emphasized 
that ESD of q\q2 is always accompanied with ESB of 
r\r2' However, this is not always true. To clarify this, we 
examine the condition to obtain ESD of the qubits and 
the companying ESB of the reservoirs. From Eqs. ([7]) 
and dHI it is obvious that the condition is Qq^q^(t) < 
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FIG. 5: (Color online) Entanglement evolution for the Ohmic 
spectral density. The two sets of parameters (77, A) = 
(0.1,5cl;o) and (0.3, lOcJo) have been considered for compar- 
ison. The corresponding entanglement evolutions are given in 
(b) and (c), respectively. In both cases a — 0.55. 



and Qrir2{t') > at any t and t' ^ which means 



\b{t')f<\a\/VT^W<l-\mf 



(14) 



When the bound states is absent, |&(oo)| = 0, the con- 
dition ([H]) can be satisfied when a < 1/a/2. So one 
can always expect ESD of the qubits and the company- 
ing ESB of the reservoirs in the region \a\ < 1/a/2, as 
shown in Fig. [2] and Refs. [19|, l32|, [33| . However, when 
the bound states are available, the situation changes. In 
particular, when \b{t)\ > ^ in the full range of time evo- 
lution, no region of a can make the condition (p!4|) to 
be satisfied anymore. For clarification, we present three 
typical behaviors of the entanglement distribution in Fig. 
m In all these cases the bound states are available. Fig. 
m^a) shows the situation where the entanglement is stably 
distributed among all of the bipartite subsystems. Fig. 
m^b) indicates that the entanglement of rir2 shows ESB 
and revival, while the entanglement of qiq2 does not ex- 
hibit ESD. Fig. m^c) shows another example that while 
the entanglement of qiq2 has ESD and revival [8[, the 
entanglement of rir2 does not show ESB but remains to 
be zero. Both Fig. |ll^b,c) reveal that ESD in qiq2 has no 
direct relationship with ESB in rir2. 



IV. DISCUSSION AND SUMMARY 

The above discussion is not dependent of the explicit 
spectral density of the individual reservoir. To confirm 



this, we consider the reservoir in free space. The spec- 
tral density has the Ohmic form J{uj) = 7^a;exp(— cj/A), 
which can be obtained from the free-space dispersion re- 
lation uj = ck. One can verify that the condition for the 
formation of bound states is: ujq — rjA < [18]. In Fig. 
[3 we plot the results in this situation. The previous re- 
sults can be recovered when the bound states are absent 
[l9|. On the contrary, when the bound states are avail- 
able, a stable entanglement is established among all the 
bipartite partitions. Therefore, we argue that the stable 
entanglement distribution resulted from the bound states 
is a general phenomenon in open quantum system when 
the non-Mar kovian effect is taken into account. 

In summary, we have studied the entanglement dis- 
tribution among all the bipartite subsystems of two 
qubits embedded into two independent amplitude damp- 
ing reservoirs. It is found that the entanglement can be 
stably distributed in all the bipartite subsystems, which 
is much different no matter to the Markovian approx- 
imate result il9il or to the decoherenceless double J-C 
model result [21|, and an identity about the entangle- 
ment in all subsystems is always satisfied. This identity 
is shown to be independent of any detail of the reservoirs 
and their coupling to the qubit, which affect only the ex- 
plicit time evolution behavior and the final distribution. 
The result is significant to the study of the physical na- 
ture of entanglement under decoherence. It implies an 
active way to protect entanglement from decoherence by 
modifying the properties of the reservoir via the poten- 
tial usage of the newlyemerged technique, i.e. quantum 
reservoir engineering [36|, [33 • 
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